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A new cosmological model based on the de Sitter gauge theory (dSGT) is studied in this paper. 
By some transformations, we find, in the dust universe, the cosmological equations of dSGT could 
form an autonomous system. We conduct dynamics analysis to this system, and find 9 critical 
points, among which there exist one positive attractor and one negative attractor. The positive 
attractor shows us that our universe will enter a exponential expansion phase in the end, which is 
similar to the conclusion of ACDM. We also carry out some numerical calculations, which confirms 
the conclusion of dynamics analysis. Finally, we fit the model parameter and initial values to the 
Union 2 SNIa dataset, present the confidence contour of parameters and obtain the best-fit values 
of parameters of dSGT. 

^ ■ 

o 

(N 

> 
O 

■ At the end of the last century, the astronomical observations of high redshift type la supernovae (SNIa) indicated 
OO ■ that our universe is not only expanding, but also accelerating, which conflicts with our deepest intuition of gravity. 

With some other observations, such as cosmic microwave Background radiation (CMBR), baryon acoustic oscillations 
' Q ' (BAG) and large-scale structure (LSS), physicists proposed a new standard cosmology model, ACDM, which introduces 
O^' the cosmological constant back again. Although this unknown energy component accounts for 73% of the energy 
density of the universe, the measured value is too small to be explained by any current fundamental theories. [l|-j3| 
bJQ, If one tries to solve this trouble phenomenologically by setting the cosmological constant to a particular value, the 
so-called fine-tuning problem would be brought up, which is considered as a basic problem almost any cosmological 
model would encounter. A good model should restrict the fine-tuning as much as possible. In order to alleviate this 
problem, various alternative theories have been proposed and developed these years, such as dynamical dark energy, 
modified gravity theories and even inhomogeneous universes. Recently, a new attempt, called torsion cosmology, has 
. attracted researchers' attention, which introduces dynamical torsion to mimic the contribution of the cosmological 
OO ■ constant. It seems more natural to use a pure geometric quantity to account for the cosmic acceleration than to 
introduce an exotic energy component. 

Torsion cosmology could be traced back to the 1970s, and the early work mainly focused on issues of early universe, 
such as avoiding singularity and the origin of inflation. In some recent work, researchers attempted to extend the 
investigation to the current evolution and found it might account for the cosmic acceleration. Among these models, 
Poincare gauge theory (POT) cosmology is the one that has been investigated most widely. This model is based on 
POT, which is inspired by the Einstein special relativity and the localization of global Poincare symmetry Goenner 
et al. made a comprehensive survey of torsion cosmology and developed the equations for all the PGT cases. [E^ Based 
on Goenner 's work, Nester and his collaborators Q found that the dynamical scalar torsion could be a possible reason 

■ for the accelerating expansion. Li et al.[7] extended the investigation to the late time evolution, which shows us the 
fate of our universe. 

Besides PGT cosmology, there is another torsion cosmology, de Sitter gauge theory (dSGT) cosmology, which can 
also be a possible explanation to the accelerating expansion. This cosmological model is based on the de Sitter 
gauge theory, in which gravity is introduced as a gauge field from de Sitter invariant special relativity (dSSR), via 
the localization of de Sitter symmetry. [8| dSSR is a special relativity theory of the de Sitter space rather than the 
conventional Minkowski spacetime, which is another maximally symmetric spacetime with an uniform scalar curvature 
1/R. And the full symmetry group of this space is de Sitter group, which unifies the Lorentz group and the translation 
group, putting the spacetime symmetry in an alternatively interesting way. But in the limit of ^ oo, the de Sitter 
group could also degenerate to the Poincare group. To localize such a global symmetry, de Sitter symmetry, requires us 
to introduce certain gauge potentials which are found to represent the gravitational interaction. The gauge potential for 
de Sitter gauge theory is the de Sitter connecion, which combines Lorentz connection and orthonormal tetrad, valued 
in 50 (1,4) algebra. The gravitational action of dSGT takes the form of Yang-Mills gauge theory. Via variation of the 
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action with repect to the the Lorentz connection and orthonormal tetrad, one could attain the Einstein-hke equations 
and gauge- Hke equations, respectively. These equations comprise a set of complicated non-linear equations, which are 
difficult to tackle. Nevertheless, if we apply them to the homogeneous and isotropic universe, these equations would 
be much more simpler and tractable. Based on these equations, one could construct an alternative cosmological model 
with torsion. Analogous to PGT, dSGT has also been applied to the cosmology recently to explain the accelerating 
expansion. [9] 

Our main motivation of this paper is to investigate (i)whether the cosmological model based on de Sitter gauge 
theory could explain the cosmic acceleration; (ii) where we are going, i.e., what is the fate of our universe; (iii) the 
constraints of the parameters of model imposed by means of the comparison of observational data. By some analytical 
and numerical calculations, we found that, with a wide range of initial values, this model could account for the current 
status of the universe, an accelerating expanding, and the universe would enter an exponential expansion phase in 
the end. 

This paper is organized as follows: First, we summarize the de Sitter gauge theory briefly in SecUll and then show 
the cosmological model based on de Sitter gauge theory in Sec. IIIIl Second, we rewrite these dynamical equations as 
an autonomous system and do some dynamical analysis and numerical discussions on this system in the Sec.lTVlandlVl 
Next in the I Vlt h section, we compare the cosmological solutions to the SNIa data and constrain the parameters. Last 
of all, we discuss and summarize the implications of our findings in Section [Vll[ 



II. DE SITTER GAUGE THEORY OF GRAVITATION 



In dSGT, the de Sitter connection is introduced as the gauge potential, which takes the form as 
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where B^^^ = r]^^ B^q^^ B^\ = r]^^ B^q^ and r]^^ = diag(l, —1,-1,— 1,-1), which combines the Lorentz connec- 
tion and the orthonormal tetrad ^ . The associated field strength is the curvature of this connection, which is defined 
as 
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where e\^^ = e^^ei 
of Lorentz-connection, 



e^e6/i,ea/i = Vab^^^, R is the de Sitter radius, and F^^^^ and T^^ are the curvature and torsion 
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which also satisfy the respective Bianchi identities. 

The gauge-like action of gravitational fields in dSGT takes the form, Q 
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Here, e = det(eJJ), ^ is a dimensionless constant describing the self-interaction of the gauge field, x is a dimensional 
coupling constant related to g and R, and F = —\F^^^e^^ is the scalar curvature of the Cartan connection. In order 
to be consistent with Einstein-Cartan theory, we take x — l/(167rG) and hg~^ = 3xA~^, where A = 3/R^. 

Assuming that the matter is minimally coupled to gravitational fields, the total action of dSGT could be written 
as: 



St = Sg ^ Sm, 



(7) 



In this paper, the Greek indices, fi,u, are 4D coordinate indices, whereas the capital Latin indices A, B, C, and the lowercase Latin 
indices, a,b,..., denote 5D and 4D orthonormal tetrad indices, respectively. 



3 



where Sm denotes the action of matter, namely the gravitational source. Now we can obtain the field equations via 
variational principle with respect to e^^B^^^^ 

-WttGx {elT,^^T\^ - ^e^^T.^^T^^) = i^GT^^ (8) 
V.F^r - {y^.ej'^ + y\.<^at + 2r[:'e,],) = 167rGi?-25M,^, (9) 



where 
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^Ma g ' ^Mab - 2^5Bf' ^^^^ 

represent the effective energy- momentum density and spin density of the source, respectively, and 

Y\,:=liT\^+T^\+T,\), (11) 

is the contorsion. It is worth noticing that the Nabla operator in Eq. (|8]) and ([9]) is the covariant derivative compatible 
with Christoffel symbols {^^} for coordinate indices, and Lorentz connection B^^ for orthonormal tetrad indices. 
Readers can be referred to Ref.[9] for more details on dSGT. 



III. THE COSMOLOGICAL EVOLUTION EQUATIONS 



Since current observations favor a homogeneous, isotropic universe, we here work on a Robertson- Walker (RW) 
cosmological metric 
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For Robertson- Walker metric, the nonvanishing torsion tensor components are of the form ^, 

T%{t)=T+{t) 5], T;,it)=T.it) e^fe, 



(12) 



(13) 



where denotes the vector piece of torsion, namely, in components, the trace of the torsion, and T_ indicates the 
axial-vector piece of torsion, which corresponds in components to the totally antisymmetric part of torsion. T+ and 
T_ are both functions of time t, and their subscripts, + and -, denote the even and odd parities, respectively. 
The nonvanishing torsion 2-forms in this case are 
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where = dt, = -^^r^ = a{t)rdO and = a{t)r smOdcj). 



^ Here, the Latin indices i, j, k..., are 3D orthonormal tetrad indices with range 1, 2, 3. 
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According to the RW metric Eq.([T2j) and the torsion Eq. ([13]), the field equations could be reduced to 
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where Eqs. ffTS]) and (p!6]) are the (t,t) and (r, r) component of Einstein-like equations, respectively; and Eqs. (pT|) and 
(p!8|) are 2 independent Yang-like equations, which is derived from the (r, ^, (/)) and (t, r, r) components of Lorentz 
connection. The spin density of present time is generally thought be very small which could be neglected. Therefore, 
we here assumed the spin density is zero. 

The Bianchi identities ensure that the energy momentum tensor is conserved, which leads to the continuity equation: 
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Equation ([T9|) can also be derived from Eqs. (p!5]) - (p!8|) . which means only four of Eqs. (p!5]) - ([19]) are independent. 
With the equation of state (EoS) of matter content, these four equations comprise a complete system of equations for 
five variables, a(t), T+(t), T_(t), p(t) and p(t). By some algebra and differential calculations, we could simplify these 
5 equations as: 
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where H = d/a is the Hubble parameter. 
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Critical Points 


Eigenvalues 
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TABLE I: The critical points and their corresponding eigenvalues. The point 9 is not physically acceptable, for its 

negative energy density. 



IV. AUTONOMOUS SYSTEM 



If we rescale the variables and parameters as 



loH; k ilk; R R/Iq; 
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where Iq = I/Hq is the Hubble radius in natural units, these variables and parameters would be dimensionless. Under 
this transformation, Eqs. (|2Q|) - (|24|) remain unchanged expect for the terms including AirGp/S and 47rGp/3, which 
change into p and p respectively. The contribution of radiation and spatial curvature in current universe are so small 
that it could be neglected, so we here just consider the dust universe with spatial flatness, whose EoS is equal to 
zero. By some further calculations, these equations could be transformed to a set of six one-order ordinary derivative 
equations, which forms a six-dimensional autonomous system, as follows. 
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For such an autonomous system, we can use the dynamics analysis to investigate its qualitative properties. Critical 
points are some exact constant solutions in the autonomous system, which indicate the asymptotic behavior of 
evolution. For example, some solutions, such as heteroclinic orbits, connect two different critical points, and some 
others, such as homoclinic orbits, are a closed loop starting from and coming back to the same critical point. In the 
dynamics analysis of cosmology, the heteroclinic orbit is more interesting. [lo| Thus, critical points could be treated as 
the basic tool in dynamics analysis, form which one could know the qualitative properties of the autonomous system. 
By some algebra calculation, we find all the 9 critical points (i^c, Qc-, ^+C7 T-ci Pc) of this system, as shown 
in Table 1. Furthermore, we analyze the stabilities of these critical points by means of the first-order perturbations. 
Substituting these linear perturbations into these dynamcial equations, we would obtain the perturbation equations 
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TABLE II: The stability properties of critical points 
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FIG. 1: The {H^T^^p) section of the phase diagram with R=4/3. The heteroclinic orbits connect the critical point 

(i) and (ii). 



around the critical points, i.e. 

5x = Ax, A=^\^=^^, (32) 

where x means the six variables of this autonomous system and / denotes the corresponding vector function on the 
right-hand side of Eqs. (|26|) - (|3T]) . Using the coefficient matrix A's eigenvalue, we could analyze the stabilities of these 
critical points. And the classification of these critical points is shown in Table HIl Among these critical points, there 
are only one positive attractor, i.e. point (i), whose eigenvalues are all negative, and only one negative attractor, i.e. 
point (ii), whose eigenvalues are all positive. The negative attractor works as a source, from which the phase orbits 
start off, whereas the positive attractor works as a sink, which the orbits finally approach. And it is the heteroclinic 
line that connects the positive attractor and the negative attractor, as shown in Fig. 1. Positive attractors are stable 
exact solutions, describing the infinite future behavior of evolution, while the unstable negative attractors depict the 
stories of infinite past. Therefore the positive attractor, point (i), here shows us the picture of late time universe, 
where all quantities tend to zero, except the Hubble parameter which approaches a finite value. At that time, the 
whole universe is entering an exponential expansion phase, just like the ACDM model. 



V. NUMERICAL DEMONSTRATION 



In order to confirm these qualitative results derived from dynamics analysis and know better about the global 
properties of this model, we explore the autonomous system by numerical methods. To solve the Eqs. (|26]) - (|3T]) 
numerically, we choose some generic initial conditions and parameters, as shown in Table IIIII First, we vary initial 
conditions (Pq, Qo, ^-o, ^+0, Po) with a fixed de Sitter radius, and the results are shown in Fig. 2(a). Then we change 
the de Sitter radius, and show the results in Fig. 2(b). Because of the rescaling Eqs. ([25]) . the current Hubble parameter 
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TABLE III: The values of initial conditions and parameters for the evolution curves in Fig. [2l 
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FIG. 2: The evolution of Hubble parameter H with respect to some initial values and parameter choice 
{R^Hq^Pq^Qq^Tj^q^ T_o^Po). According to the transformations (j25j), the unit of time here is the Hubble Time. In 
Fig 2(a), we fixed R and changed T±, while in Fig 2(b), we changed R. 




— (0.4,1,0,0,-1.5,0,0.8) 

— (0.6,1,0,0,-1,0,1) 

— (1.5,1,0,0,0,0,0.5) 



here must be 1. From these numerical results, it is easy to find that the Hubble paramter of all the solutions approaches 
a particular finite value in the infinite future, whatever the initial conditions are, and this value only depends on the 
de Sitter radius R. Such results demonstrate the dynamics analysis we have done in the former section. We could also 
find that this positive attractor covers a wide range of initial conditions, and therefore the troublesome fine-tuning 
problem has been alleviated. In comparison with the result of PGT, we find the cosmology based on de Sitter gauge 
theory is quite different from the Poincare gauge theory, where the expansion will asymptotically come to a halt. It 
is the existence of the de Sitter radius that makes such a discrepancy. If we let ^ oo, the de Sitter gauge theory 
would degenerate to the PGT. 



VI. SUPERNOVAE DATA FITTING 



A basic approach to testing a cosmological model is the supernova fitting through its description of the expansion 
history of the universe. In this section we fit the initial conditions and model parameters to current type la supernovae 
data. And the maximum likelihood technique is used here, which could determine the best fit values of parameters 
and initial conditions and the goodness of this model. The supernova data are comprised of the distance modulus 
fiobs, which is equal to the difference between apparent magnitude rrii and absolute magnitude M^, and redshifts Zi of 
supernovae with their corresponding errors (Ji. Note that the error here are assumed to be normally distributed and 
independent. 

The theoretical distance modulus is related to the luminosity distance dL hy 

„tu{zi) = 51ogio(^)+25 

= 51ogioi?L(^i)-51ogio(^)+25 

= 51ogio/)L(^i)-51ogio/i + 42.38, (33) 
where the Dl{z) is the dimensionless 'Hubble-constant free' luminosity defined by Dl{z) = HodL{z)/c. 
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For a spatially fiat cosmological model, which we consider here, the luminosity distance could be expressed in terms 
of Hubble parameter H{z)^ as follows, 

Dl{z) = (1 + z) Td/-— i (34) 

Jo H{z';ai,...an) 

where the Hubble parameter H{z^; ai^ ...an) here is the dimensionless Hubble parameter under the rescaling transfor- 
mation Eq ([25]) . 

As we know, due to the normal distribution of errors, we could use the parameter as the maximum likelyhood 
estimator to determine the best fit values of parameters and initial conditions (i?, Pq, Qo, ^+o, ^-o, Po) of the model. 
The here for the SNIa data is 

X W = 2^ , 

i * 

^ J2 ll^obsi^i) -^^ogioDL{zi;0) - fip]^ ^^^^ 

i * 

where jhq = — Slog^o h-\- 42.38 , denotes all the parameters and initial conditions, and are the statistical errors of 
SNIa. If we want to include the systematic errors, which are comparable to the statistical errors and should be taken 
into account seriously, we could resort to the covariance matrix Csn^ and the Eq. (|35]) turn out to take the form 

N 

= i^^obs{Zi) - fith{Zi)] {CsN)ij[f^obs{Zj) - fJ^thiZj)], (36) 

N 

= Y il^obsizi) - 5\og^QDL{zi]0) - /io] {Cg^)ij[/j.obs{zj) - 5\og^QDL{zi;0) - /io]. (37) 

The parameter /io here is a nuisance parameter, whose contribution we are not interested in. So we marginalize over 
this parameter, /io, thus obtaining a new 



x\e) = A{0)-^ + ln[-\, (38) 



where 



N 



MO) = Yl [^'obs{z^; 0) - 5 logio DL{z^; 0)] {Cs]:,)^J[^iobs{zJ;0) - 5 logio DL{zj;0)i (39) 

N 

B{9) = Y.{Csl,)ij[^iobs{zf,0)-^\og,^DL{zf,9)l (40) 

3 

N 

C{9) = Y, = {Cslj)ii. (41) 



Now we try to constrain our model parameter and initial values by this maximum likelihood estimator. The 
dataset we use here is the "Union2" SNIa dataset (N=557), the most comprehensive one to date, which combines all 
the previous SNIa dataset in a homogeous manner. 

By minimizing the we find the best fit parameters of dSGT model, as shown in Table. HVl Based on the current 
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T+o 
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1.1005 














0.0431 535.3384 



TABLE IV: The best-fit initial data and parameters 



observations, the present density of torsion in our universe is very small, so it is reasonable to assume that the initial 
values of all torsions and their first order derivative are zero at z = 0. But their second order derivatives does not 
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FIG. 3: The distribution with respect to R and p, compared to the ACDM, the plane = 536.634. Here we 
assume that all the torsions and their first order derivatives vanish at present time. 



0.4 
0.35 
0.3 
0.25 
Q_ 0.2 
0.15 
0.1 
0.05 


1 1.02 1.04 1.06 1.08 1.1 1.12 1.14 1.16 1.18 

R 

FIG. 4: The 68.3%, 95.4% and 99.7% confidence contours of dSGT with respect to R and p, using the Union 2 
dataset. Here we also assume that the current-time values for all the torsions and their first order derivatives are 

zero. The yellow point is the best-fit point. 

vanish yet, which would have a significant impact on the history and future of the evolution of our universe. In this 
case, the number of parameter and initial value is reduced, and the rest parameters and initial values are just R and 
po- It is easy to find the best-fit of these 2 parameters, which are shown in Table HVl And the minimal is 535.3384, 
whereas the value for ACDM is 536.634, with Vtrn = 0.27, I^a = 0.73. In Fig. [3] we show the distribution with 
respect to R and p compared to ACDM model, where the plane x^ = 536.635 corresponds to the value of ACDM. 
Furthermore, we plot the contours of some particular confidence levels, as shown in Fig. HJ From these figures, we 
could find that the evolution of our universe is insensitive to the initial value, which alleviate the fine-tuning problem. 



546 




10 



VII. SUMMARY AND CONCLUSION 



The astronomical observations imply that our universe is accelerating to a de Sitter spacetime. This gives us a 
strong motive to consider the cosmic evolution based on the de Sitter gauge theory instead of other gravity theories. 
The localization of de Sitter symmetry requires us to introduce curvature and torsion. So in de Sitter gauge theory, 
the torsion is an indispensable quantity, by which people tried to include the effect of spin density in gravity theory 
at first. But now this essential quantity might account for the acceleration of our universe, if we apply dSGT to 
cosmology. 

We found the cosmo logical equations for dust universe in dSGT could form an autonomous system by some transfor- 
mations, where the evolution of the universe is described in terms of the orbits in phase space. Therefore, by dynamics 
analysis to the dSGT cosmology, one could study the qualitative properties of this phase space. We found all 9 critical 
points, as shown in Table HI We also analyzed the stabilities of these critical points, and found among these critical 
points there is only one positive attractor, which is stable. The positive attractor alleviates the fine-tuning problem 
and implies that the universe will expand exponentially in the end, whereas all other physical quantities will turn 
out to vanish. In this sense, dSGT cosmology looks more like the ACDM,than PGT cosmology. And we conducted 
some concrete numerical calculations of this the destiny of our model of the universe, which confirms conclusions from 
dynamics analysis. 

Finally, in order to find the best-fit values and constraints of model parameters and initial conditions, we fitted 
them to the Union 2 SNIa dataset. The maximum likelihood estimator here we used is the estimate. By minimizing 
the we found the best-fit parameters R = 1.135, p = 0.274 and the corresponding = 535.3384, while the value 
for ACDM is 536.634, with Qrn = 0.27, Qa = 0.73. Note that we here set all the initial values of torsions and their 
first-order derivatives to zero at t = to, since the contribution of torsion to the current universe is almost negligible. 
We also plotted the confidence contour Fig.HJwith respect to R and p, from which it is easy to see that the fine-tuning 
problem is alleviated and the evolution is not so sensitive to the initial values and model parameters. 

If we want to go deeper into cosmology based on de Sitter gauge theory, there are a lot of work need to be done. 
We should fit this model to some other observations, like BAO and LSS etc, to constrain the parameters better. We 
also could study the perturbations in the early universe, and compare the results to CMBR data. These issues will 
considered in some upcoming papers. 



Acknowledgments 

This work is supported by SRFDP under Grant No. 200931271104 and Shanghai Natural Science Foundation, 
China, under Grant No. 10ZR1422000. 



[1] P. J. E. Peebles and B. Ratra, The cosmological constant and dark energy, Rev. Mod. Phys. 75 (2003) 559-606, 
I astro-ph /0207347 . 

[2] T. Padmanabhan, Cosmological constant: The weight of the vacuum, Phys. Kept. 380 (2003) 235-320, hi5^th702T2290l ; 

X.-Z. Li and J.-G. Hao, 0(N) phantom, a way to implement w < -1 Phys. Rev. D69 (2004) 107303, hep-th/0303093 . 
[3] J.-G. H ao and X. -Z. Li, Phantom Cosmic Dynamics: Tracking Attractor and Cosmic Doomsday, Phys. Rev. D70 (2004) 

043529, [astro-th/0309746 ; 

C.-J. Feng and X.-Z. Li, Cardassian Universe Constrained by Latest Observations Phys. Lett. B692 (2010) 152-156, 
iarXiv:0912.4793^. 

[4] F. W. Hehl, Four lectures on Poincare gauge field theory in Proc. of the 6th Course of the School of Cosmology and 
Cravitation on Spin, Torsion, Rotation, and Supergravity, eds. P. G. Bergmann and V. De Sabbata (New York: Plenum) 
p5 

M. Blagojevic, Cravitation and Cauge Symmetries, loP Publishing, Bristol, 2002. 
[5] H. Goenner and F. Miiller-Hoissen Spatially homogeneous and isotropic spaces in theories of gravitation with torsion. Class. 

Quantum Crav. 1 (1984) 651-672. 
[6] H. Chen, F.-H. Ho and J. M. Nester, C.-H. Wang, and H.-J. Yo, Cosmological dynamics with propagating Lorentz connection 

modes of spin zero, JCAP 0910 (2009) 027, arXiv:0908.3323j 

P. Baekler, F. W. Hehl, and J. M. Nester, Poincare gauge theory of gravity: Friedman cosmology with even and odd parity 
modes. Analytic part, Phys. Rev. D83 (2011) 024001, arXiv:1009.5112j 

F.-H. Ho and J. M. Nester, Poincare Cauge Theory With Coupled Even And Odd Parity Dynamic Spin-0 Modes: Dynamic 
Equations For Isotropic Bianchi Cosmologies, arXiv: 1 106.071 1; 

F.-H. Ho and J. M. Nester, Poincare gauge theory with even and odd parity dynamic connection modes: isotropic Bianchi 
cosmological models, arXiv:1105.5001 . 



11 



[7] X.-Z. Li, C.-B. Sun, and R Xi, Torsion cosmolgical dynamics, Phys. Rev. D79 (2009) 027301, [arXiviOgOS.SOSSj 

X.-Z. Li, C.-B. Sun, and P. Xi, Statefinder diagnostic in a torsion cosmology, JCAP 0904 (2009) 015, arXiv:0903.4724"; 
X.-C. Ao, X.-Z. Li, and P. Xi Analytical approach of late-time evolution in a torsion cosmology Phys. Lett. B694 (2010) 
186-190, [arXiv:101.4117]. 

[8] H.-Y. Guo, C.-G. Huang, Z. Xu, and B. Zhou, On Beltrami model of de Sitter spacetime, Mod. Phys. Lett. A19 (2004) 
1701-1710, hep-th/0311156 ; 

H.-Y. Guo, C.-G. Huang, Z. Xu, and B. Zhou, On special relativity with cosmological constant, Phys. Lett.. A331 (2004) 
l-7,[hep-th/043171]; 

H.-Y. Guo, C.-G. Huang, Z. Xu, and B. Zhou, Three kinds of special relativity via inverse Wick rotation, Chin. Phys. Lett. 
22 (2005) 2477-2480, hep-th/0508094 ; 

H.-Y. Guo, C.-G. Huang, B. Zhou, Temperature at horizon in de Sitter spacetime, Europhys. Lett. 72 (2005) 1045-1051. 
[9] C.-G. Huang, H.-Q. Zhang, and H.-Y. Guo, Cosmological Solutions with Torsion in a Model of de Sitter Cauge Theory of 
Cravity JCAP 0810 (2008) 010, [arXiv:0801.0905_. 
[10] X.-Z. Li, Y.-B. Zhao, and C.-B. Sun, Heteroclinic orbit and tracking attractor in cosmological model with a double expo- 
nential potential. Class. Quant. Crav. 22 (2005) 3759-3766, [astro-ph/050819]; 

J.-G. Hao and X.-Z. Li, An attractor solution of phantom field, Phys. Rev. D67 (2003) 107303, |gr- qc/0302100] . 



